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Abstract 


Let  the  positive  definite  matrix  A have  a Cholesky  factorization 
A - R k.  For  a given  vector  x suppose  that  A ■ A - xxT  has  a 
Oiolesky  factorization  A - RTR.  This  paper  considers  an  algorithm 
for  confuting  R from  R and  x and  an  extension  for  removing 
a row  from  the  QR  factorization  of  a regression  problem.  It  is 
shown  that  the  algorithm  is  stable  in  the  presence  of  maiding 
errors.  However,  it  is  also  shown  that  the  matrix  R can  be  a 
very  ill-conditioned  fine t ion  of  R and  x. 
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The  Effects  of  Rounding  Error  on  an  Algorithm 
for  Downdating  a Qiolesky  Factorization 


G.  W.  Stewart 


1.  Introduction 


',rxrs$' 

V *v 

/ 

fins 

. 'DC 

niv  ■ 

JIX. 

sr 

□ 

BY 

DISTniB 

, coa 

Oisl.  •• 

'1...  1, 

i.'iCIAL 

A 

Let  A be  a positive  definite  matrix  of  order  p.  Then  A can 


be  factored  in  the  form 


A = 


where  R is  upper  triangular.  This  "Qiolesky  factorization"  of  A is 
unique  up  to  the  signs  of  the  rows  of  R (e.g.  see  [61 

In  this  paper  we  shall  be  concerned  with  the  following  problem. 

Given  a p-vector  x and  the  matrix  R find  the  Qiolesky  factorization  of 
the  matrix 

(1.1)  A - A - xxT  , 

where  it  is  assumed  that  x is  such  that  A is  positive  definite.  We 
shall  refer  to  this  problem  as  the  downdating  problem. 

An  important  application  of  downdating  is  the  removal  of  an  observa- 
tion from  a linear  regression  problem  that  is  being  solved  by  means  of  the 
QR  factorization.  Specifically,  consider  the  problem  of  minimizing 

p2  - HXp-yll2  , 

where  X is  an  n * p matrix  of  rank  p and  ||*||  denotes  the  usual 

2 T 

Euclidean  vector  norm  defined  by  i|x||  - x x.  It  is  well  known  that  X 
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can  be  factored  in  the  form 

X - QR  , 

where  R is  upper  triangular  and  Q has  orthonormal  colunns,  i.e. 

QTQ  * I.  If  z is  defined  by 

z ■ QTr  , 

then  the  solution  of  the  regression  problem  is  given  by 

p - R Xz 

and  the  residual  sum  of  squares  by 

p2  - llyll2-ll  zll 2 . 

When  n is  large,  it  may  be  inq>ossible  to  retain  the  elements  of 
the  n * p matrix  Q in  the  main  memory  of  the  computer  performing  the 
calculation.  In  this  case  one  may  compute  R,  z,  and  p without  expli- 
citly forming  Q [3,5].  Although  this  suffices  for  the  computation  of 
p,  one  is  left  with  the  problem  of  performing  a variety  of  statistical 
computations  when  one  knows  only  R,  z,  and  p.  (It  is  interesting 

to  note  that  aficiandos  of  the  normal  equations  have  the  same  problem; 

T 

they  cannot  retain  X in  main  memory  and  must  work  instead  with  X X, 
XTy,  and  ||y||2  [1].) 

One  frequently  occurring  requirement  is  to  remove  an  observation 

T 

from  the  regression,  that  is  to  remove  a row  x from  X and  the  corre 
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sponding  component  r\  from  y.  Without  loss  of  generality  we  may 
T 

suppose  that  x is  the  last  row  of  X,  so  that  X can  be  written 
in  the  form 


It  follows  that 

~T~  T T 

(1.2)  X*X  - X*X  - xx  . 

Now 

XTX  * RTQTQR  - RTR  , 

which  shows  that  the  triangular  part  of  the  QR  factorization  of  X is 

T ^ 

the  Cholesky  factor  of  X X.  Likewise  R,  the  triangular  part  of  the  QR 
factorization  of  X,  is  the  Cholesky  factor  of  X!X.  Coirparing  (1.1)  and 

(1.2) ,  we  see  that  the  problem  is  one  of  downdating  the  Cholesky  factori- 

T 

zation  of  X X.  There  is,  of  course,  more  to  it  than  this,  for  we  must 

~ ~2 
also  compute  the  downdated  vector  z and  residual  sun  of  squares  p . 

In  this  paper  we  shall  give  a rounding  error  analysis  of  an  algorithm 
for  computing  R,  z,  and  p.  Our  conclusions  are  that  the  algorithm  is 
remarkably  stable;  however,  this  stability  does  not  guarantee  that  the 
results  are  accurate,  for  the  downdating  problem  can  be  quite  ill  condi- 
tioned. We  begin  with  a discussion  of  this  ill  conditioning,  before  going 
on  to  a description  of  the  algorithm  and  the  subsequent  error  analysis. 
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2.  The  condition  of  the  downdating  problem 

Let  A,  A,  and  x be  as  in  the  previous  section  with  A having 

T 

a Cholesky  factorization  R R.  Let  A have  the  Cholesky  factorization 
R R.  In  applications  we  will  of  course  not  know  A and  B.  Rather  we 

will  be  given  R and  x and  be  required  to  compute  ft.  Consequently, 

we  are  interested  in  assessing  the  effects  of  perturbations  in  R and 
x on  R. 

We  shall  first  consider  a perturbation  E in  R.  Assune  that  the 
T T 

matrix  (R+E)  (R+E)  - xx  has  a Cholesky  factor  R.  We  wish  to  assess 
the  size  ||R-R||,  where  here  INI  denotes  the  spectral  matrix  norm 
[6,7],  We  begin  by  comparing  the  singular  values  [6]  of  R and  R, 
which  we  denote  by  o^  > - •••  i 8 and  ^ - • • • - Op.  Now 

RTR  - (R+E) T (R+E)  - xxT 

» RTR  - xxT  ♦ RTE  + ETR  + ETE 

= R^R  + RTE  + ETR  + ETE  . 

222  t «2  • 2 -2 

Since  are  the  eigenvalues  of  R R and  likewise  of ,o ,0'; 

I P A L p 

-T- 

are  the  eigenvalues  of  R R,  it  follows  from  the  classical  perturbation 
theory  for  eigenvalues  of  symnetric  matrices  [6,7]  that  for  i ■ l,2,...,p 

|o^?|  * ||RTE+ETR+ETE||  < + e2  , 

where  Oj  - ||R||  and  e - ||E||.  In  particular 


Tz 

CTi 


\*t)  »i  ' 


and  it  follows  from  the  inequality  |1  - /I  + x | < |x|  that 

(2.1) 


2o1e+e^  _ 20,6+e^ 

a, , ~ CTi  5 CTi  + —71 


~ZT 

ai 


Now 

||R  - R||  > max|5i  - ct-J ; 

hence  (2.1)  has  the  disturbing  implication  that  ||R-R||  can  be  as  large  as 

2 ' 1 

(2cije+e  )/5p.  In  particular  if  5pS/2o^e+e  , we  cannot  guarantee  that 

5p  and  Op  agree  in  any  significant  figures. 

Casting  the  results  in  terms  of  relative  errors  (e.g.  rounding  errors) 

may  make  this  clearer.  Suppose  that  the  original  matrix  R has  nonzero 

elements  all  of  about  the  same  size,  and  these  elements  are  perturbed  by 

a relative  error  of  order  eM-  Then  in  the  above,  e » e^lRII  “ Enfi»  so 

that  if 

(2.2) 


*p“  M 


3p  may  be  obliterated  by  the  error  in  R.  The  square  root  has  the  impli- 
cation that  in  downdating  one  cannot  tolerate  a spread  of  singular  values 
of  half  the  computational  precision  without  losing  all  precision  in  the 


llest  singular  value. 


- 6 - 


Perturbations  in  x have  much  the  same  effect.  If 
RTR  = RTR  - (x+f)(x+f)T, 

» 

then  a repetition  of  the  above  argument  shows  that 

I,  | < ziixinifii-iifii2 . 

07 

1 

.1 

It  should  be  observed  that  the  dependence  of  the  bounds  on  5^ 

can  be  removed  by  the  following  argument.  The  derivation  of  (2.1)  is 

symmetric  in  o^  and  o^.  Consequently  we  may  replace  the  denominator 
2 2-2 

in  (2.1)  by  = max{o^,o^}.  But  we  also  have  the  bound  s 

Combining  the  two  bounds  gives 

(2.3)  |5i-oi|  < (2o16+e2)1/f3  . 

In  conputational  practice  it  is  unlikely  that  both  o^  and  o ^ will  be 
2 1/2 

less  than  (2o1e+e  ) , so  that  the  cube  root  in  (2.3)  is  effectively 

a square  root. 

That  the  bound  (2.1)  is  realistic  can  be  seen  by  considering  the 
scalar  case  p - 1.  This  case  actually  arises  in  practice;  for  in  the 
regression  problem  mentioned  in  §1,  X becomes  an  n-vector,  x becomes 
a component  of  X,  and  R becomes  ||X||.  Thus  the  downdating  problem 
becomes:  given  the  norm  of  a vector  find  the  new  norm  after  a component 

has  been  removed  from  the  vector.  The  results  of  this  section  have  the 
following  implications  for  downdating  norms  in  t-digit  arithmetic.  If 

t/2 

ever  a sequence  of  downdates  reduces  the  norm  by  a factor  greater  than  10  , 

the  results  can  be  expected  to  be  conpletely  spurious. 


3.  The  algorithm 
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The  algorithm  described  in  this  section  is  an  extension  of  one  that 
has  appeared  in  [2,5]. 

We  shall  use  the  notation  introduced  in  §1.  We  assume  that 

the  reader  is  familiar  with  computations  with  plane  rotations  (for  details 

see  [6]  or  [7]).  In  order  that  the  several  computational  steps  of  the 

algorithm  will  not  be  lost  in  the  derivations,  we  proceed  immediately  to 

a description  of  the  entire  algorithm. 

T T 

lve  the  system  a R = x . 

T T 

a||  > 1,  report  R R - xx  indefinite  and  stop. 

t " ■■■  "TW 

Compute  a = /l-||a|p  . 

4.  For  i = p,p-l,...,l  determine  plane  rotations  lb  in  the 
(i,p+l)  plane  such  that 


To  justify  this  algorithm,  we  first  show  that  the  condition  ||a||  < 1 

T T 

is  necessary  and  sufficient  for  R R - xx  to  be  positive  definite.  In 
fact 

(3.1)  (RTR-xxT)  =*  RT(I-aaT)R  . 

T 2 

Now  the  eigenvalues  of  I-aa  are  1 - ||a||  of  multiplicity  unity  and  1 

T t 

of  multiplicity  p-1.  It  follows  that  I-aa  , and  hence  ft  ft,  is  positive 


definite  if  and  only  if  |j  a||  ^ < 1. 
Let  Q * ui***Up-iUp-  Then 


(3.2) 


a R 

a 0 


0 R 


P o 


in  which  we  must  verify  that  p = 1 and  b = x.  Since  Q Q = I,  if  each 
side  of  (3.2)  is  multiplied  by  its  transpose,  the  result  is 


RTa  RTR 


~T~  T 
pb  R R+bb 


It  follows  immediately  thatp  = 1,  b * x,  and 


T Tn  T 
R*R  - RlR  + xx1 


But  it  is  easily  seen  from  the  form  of  the  plane  rotations  IL  that  R 
is  upper  triangular.  Hence  R is  the  downdated  Cholesky  factor. 

In  applications  to  regression  problems,  it  is  necessary  to  confute 
z and  p.  One  way  of  approaching  this  is  to  observe  that  the  Cholesky 
factor  of  (X,y)T(X,y)  is 


(3.3) 


Thus  the  new  decomposition  can  be  determined  by  removing  (x  rj)  from  (3.3), 


However  we  prefer  to  use  a different  algorithm  for  two  reasons.  First  if 
one  has  several  vectors  y,  the  algorithm  must  be  repeated  for  each  one. 
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with  considerable  diseconomies  in  time  and  storage.  Second,  the  aug- 
mented downdating  may  fail,  even  though  R by  itself  can  be  downdated, 
and  it  is  desirable  not  to  confound  these  sources  of  failure. 

In  the  description  of  our  proposed  algor ithm  below,  c^  and  s^ 
are  the  cosines  and  sines  defining  the  plane  rotations  U; . 

1.  Set  n0  = P 

2.  For  i = l,2,...,p  compute 


(3.4) 


rj . = s.z.  + c.n . . 
'i  ii  i'i-l 


3.  If  n > p stop 

. ~ /7T7 

4 . p = /p  -np 

To  show  that  these  formulas  indeed  produce  the  required  2 and  p, 
we  first  observe  that  they  are  well  defined,  since  a j1  0 implies  that 
no  c^  can  be  zero.  Now  the  two  relations  in  step  two  of  (3.4)  are  equi 
valent  to 


(3.S) 


It  follows  that 


Ci  ‘Si 


s.  c. 

1 1, 


whence 
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e ;,)  • ( 


R 1 

T 

X T] 


It  then  follows  that 


R z\7R  z 


0 p / V 0 p 


~ T~  T 'v  T«v 

R R+xx1  R*z+r)x 

t~  T T 2 

Z R-h^X  Z Z+p 


T 2 ~T~  2 2 

But  zz  + p = zz  + p +t};  hence 


R z\T/R 

0 p/  \0 


T/  R 2 


0 p 


and  2 and  p conprise  the  last  column  of  the  Cholesky  factor  of  the 
downdated  augmented  system. 

We  note  that  if  > p,  then  p is  not  large  enough  to  accomnodate 

T 

the  decrease  in  the  residual  due  to  the  deletion  of  (x  ,r)) , and  the  algo- 
rithm should  be  stopped. 
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4.  The  effects  of  rounding  error 

In  this  section  we  shall  adopt  the  conventions  and  assumptions  usual 
in  floating-point  rounding -error  analyses.  If  e is  an  arithmetic  expres- 
sion with  a specified  order  of  evaluation,  £1  (e)  will  denote  the  result 
of  evaluating  e in  floating-point  arithmetic.  We  shall  assume  that 
floating-point  multiplication  and  division  satisfy 

(^(aob)  * aob(l+e),  o * x f 

where 


Here  is  the  rounding  unit  of  the  conputer  in  question  (i.e.  eM  is 
approximately  the  largest  number  e for  which  ^(l+e)  * 1).  We  assune 
addition  and  subtraction  satisfy 

l^(a±b)  = a(l+ej)  ± b(l+e2) 

where  | e ^ | » | e 2 1 eM-  Finally  we  assume  that 

- (l+e)^a  , 

where  again  |e|  s eM«  As  is  customary,  we  ignore  problems  of  overflow  and 
underflow. 

In  order  to  simplify  our  bounds  we  shall  freely  discard  higher  order 
terms  in  eM-  For  example  (l+eM)  (1+e^  will  be  approximated  by  1 ♦ 2eM« 
Although  our  results  will  no  longer  have  the  status  of  theorems,  their 


I 


r 
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derivation  will  be  considerably  less  cluttered.  Moreover,  if  p is 
sufficiently  small,  say  pe^  < .01,  then  the  bounds  can  be  made  rigorous 
by  multiplying  by  a factor  near  unity. 

We  begin  with  the  computation  of  a.  Here,  and  in  what  follows, 
all  quantities  stand  for  their  confuted,  not  their  true,  values.  The 
solution  of  triangular  systems  has  been  analyzed  elsewhere  [6,7],  and 
we  merely  quote  the  results.  The  vector  a satisfies 

(4.1)  aT(R+F)  - xT 
where 

i 

(4.2)  |ftJl  = 0*2)|rtjUM  • 

It  follows  that  if  r^  and  f^  denote  the  j-th  columns  of  R and  F,  then 

(4.3)  IlfjH  < ^(j^Hr.lle^,  . 

2 

We  turn  now  to  the  computation  of  a.  If  we  compute  a in  the 

2 2 2 

order  1 - (8^32+...^)  we  have 

' a2  * (1+£q)  ' ai(1+ei)  * a2^+e2^  ap(1+ep) 


where  |eQ|  £ eM  and  | e i | £ (p-i+3)eM  (i>l).  Hence,  since  ||a||‘  < 1, 
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The  computation  and  application  of  plane  rotations  has  been  analyzed 
in  detail  by  Wilkinson  [6] , where  he  shows  that  there  are  exact  rotations 

A A A 

U^.U^, . . . ,Up  such  that  for  any  vector  v 


^U1"*UP-1UPV) 


U. . . .U 
1 P 


,U  v + 

-1  p 


where 


(4.4)  llgll  < 6p||vl! . 

Here  we  have  suppressed  some  second  order  terms  that  account  for  the  slow 
growth  in  a bound  on  ||U,  ...U  ,U  v||. 

1 p-1  p 

Let 


(4.5)  (J  - Ui...Up_itfp  . 

We  first  consider  the  application  of  0 to  the  vector  (aT,a)T.  From  the 
results  quoted  above 


(4.6) 

where 


Il80ll  < 6?^  ||a||2>a2) 1/2  . 


Now  from  (4.4) 


Hall2  ♦ a2  - ||a||2  ♦ (W^d-llall2*^) 
- 1 ♦ 2||a||Tj  ♦ t2  • 


Hence 
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(4.7) 
and 

(4.8) 


where 

(4.9) 


» M < 6peM 


P « (l!a||2+a2-|lg0l|2)1/2 


1 + CT0  * 


N £ *Tr  cm 


T T 

We  next  consider  the  application  of  Q to  (rj,0)  . We  have 


(4.10) 


where 


(4.11) 


<;)•  (I::;)  • 


l!gj  II » IYj  I < 6p|| Tj  ||eM  . 


Here  is  the  computed  value.  We  wish  to  find  a bound  on  lxj'Sjl 
Since  0 is  orthogonal,  we  have  from  (4.6),  (4.8)  and  (4.10)  that 


- (gj. 


„(V‘A 

VjV 


1 * rj  * !0rj 

‘ li  * Vj  * Tj  * #J  • 
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But  from  (4.1) 


T T 

» W “ f j ' 


Hence 


V ‘j  ■ °<f-i  * Ti  * gfo  ‘ aTfj  • 


Since  up  to  terms  of  order  ew,  ||r.||  - It  CK  ,5  Oil , we  have  from  (4.3), 

M J J J 

(4.7),  (4.9),  and  (4.11)  that 


W °j 


where 


(4.12)  | o.|  < [ Up-  * /J  (j+2)  ]||r.||eM  . 

To  suimarize  we  have  shown  that  there  is  an  orthogonal  matrix  Q 
such  that 


where  G and  s satisfy  (4.11)  and  (4.12).  In  other  words,  the  confuted 
downdated  Cholesky  factor  R is  very  near  the  factor  obtained  by  downdat- 
ing with  a slightly  perturbed  vector  x.  The  error  G in  R is  uninportant , 
except  as  it  may  affect  subsequent  downdates;  however,  the  results  of  §2 
show  that  the  error  s in  x may  seriously  affect  the  accuracy  of  R. 

Two  other  points.  First,  the  higher  order  term  in  (4.12)  is  due  to 

T T 

the  solution  of  the  triangular  system  a R ■ x . The  factor  (j+2)  can  be 
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removed  from  this  term  by  accunulating  inner  products  in  double  precision; 
however,  in  practice  this  is  unnecessary,  since  the  term  does  not  dominate 
its  companion  (and  this  only  in  colum  p)  until  p ■ 40,  and  it  is  not  yet 
double  when  p - 150. 

Second,  the  bounds  are  given  column  by  colimn  and  hence  are  indepen- 
dent of  colunn  scaling.  This  is  not  surprising,  since  the  confutations 
in  each  colunn  are  independent  of  one  another. 

We  turn  now  to  the  analysis  of  the  errors  involved  in  downdating  1. 

Define 

m ^2* * * " *^i- i»^i»  • • • * 

so  that  from  (3.5) 

Wj  “ i * 1,2, ...  ,p  . 

T 

However,  the  evaluation  of  ££(U  w^.j)  is  not  the  straightforward  one  implied 
by  (3.5);  rather  it  is  the  indirect  one  inf  lied  by  the  formulas  in  step  2 
of  (3.4),  which  we  now  analyze.  We  have 

(x  (lnOn.fl,  .(ln.Hln,)] 

\ ^ 2 — (ln4) 

where  |c^|  « e.^.  Thus 

zi  “ c^dn^^ClM,)'1  - •i^i.1(l*«2)(lns)(ln1)‘1  , 
and  it  follows  that 


2i  " ciVln5>  ‘ *iVi(ln6) 
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where  leJ  < 2e,,  and  |e,l  < 3ew.  Likewise 

= + cifii-l(1+e8)  * 

where  | e y | * | e g | 26^. 

These  results  show  that  as  far  as  rounding  errors  are  concerned, 
the  formulas  in  (3.4)  are  equivalent  to  the  direct  application  of  (3.5), 
with  the  exception  that  the  term  has  the  bound  3e^  instead  of 
2eM-  This  means  that  the  error  analysis  of  Wilkinson  cited  above  goes 
through  mutatis  mutandis,  with  the  result  that  the  right  hand  side  of  the 
bound  (4.4)  becomes  7p|| vi! eM-  Hence 

wp  - QTw0  * g , 

A 

where  ||g||  < 7p||w0lle^.  Since  Q is  orthogonal. 


where 

M,M  < 7p(||z|!  V)eM  . 

Thus  the  computed  1 is  very  near  the  vector  that  would  be  obtained 
by  downdating  z with  a slightly  perturbed  q . It  should  be  noted  that 
the  transformation  Q is  the  same  as  the  one  defined  by  (4.5)  in  the 
previous  analysis. 

It  goes  without  saying  that  these  bounds  are  an  extreme  over-estimate 


A 


of  the  errors  that  would  be  encountered  in  practice.  None  the  less  they 
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suffice  to  demonstrate  the  exceptional  stability  of  the  algorithm.  Any 
inaccuracies  observed  in  the  results  cannot  be  attributed  to  the  algorithm; 
they  must  instead  be  due  to  the  ill-conditioning  of  the  problem.  This 
raises  the  question:  does  the  algorithm  provide  some  way  of  detecting 

ill  conditioning?  We  shall  answer  this  question  in  the  next  section. 
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5.  The  meaning  of  ||a|| 

It  is  a consequence  of  the  results  of  §2  that  ill  conditioning 

in  the  downdating  problem  is  associated  with  small  singular  values  in 

R.  In  §3  it  was  shown  that  if  ||a||  * 1,  then  R is  singular,  i.e. 

5p  = 0.  It  is  therefore  reasonable  to  conjecture  that  values  of  a 

near  unity  will  be  associated  with  ill-conditioned  problems  and  vice 

versa.  However,  just  as  the  determinant  is  a poor  indicator  of  the 

condition  of  a matrix,  the  value  of  II all  may  be  a poor  indicator  of  the 

condition  of  the  downdating  problem.  In  this  section  we  shall  show  that 

the  value  of  ||a||  will  reliably  signal  trouble. 

We  first  show  that  the  value  of  ||a||  cannot  cry  wolf;  if  it  is 

near  inity,  then  the  problem  must  be  ill  conditioned.  It  follows  from 

T 2 

(3.1)  and  the  fact  that  the  smallest  eigenvalue  of  I - aa  is  1 - ||a|| 
that  the  smallest  eigenvalue  of  R R is 

- ||R|||(l-||a||2)  . 

Since  5p  - 


2 

It  follows  from  the  discussion  surrounding  (2.2)  that  if  1 - ||a||  ■ 0(eM) 
then  R can  be  expected  to  lose  about  half  its  accuracy. 

We  cannot  show  that  a small  value  of  ap  implies  that  ||a||  is  near 
unity.  However  we  can  show  that  if  any  singular  value  of  R is  reduced 
in  the  downdating  by  a significant  factor,  then  ||a||  must  be  near  unity. 
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We  start  by  developing  an  expression  for  ||a||.  First 


I a|| 2 - xTR-1R  Tx  - xT(RTR)  "1x 


T/-dTX.  T,-l 
x (R  R+xx  ) x 


Set 


so  that 


xTR_1  (I+R^xx^R*1)^"^  . 


b * R Tx  , 


| a|| 2 - bT(I+bbT)  '1b  . 


Since  b is  an  eigenvector  of  I ♦ bb  corresponding  to  the  eigenvalue 
1 + ||b||2,  it  follows  that 


(5.1) 


.II2  - . 

i*Bbir 


We  next  obtain  a lower  bound  on  ||b||  . Let  v^  be  the  right  singular 


vector  corresponding  to  and  let  ■ (v^  ,v^ , . . . ,Vp) . Then  if  V^x  + 0 
(5.2) 


_x  livTxii 
llbll  - HR  x||  > -i— 


But  from  the  minimax  theorems  [6,7] 

o\  < ii^rtrv.h  < iivjRTftvin  ♦ llvTxxTvin 
- o\  ♦ livjxii2  . 
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I 


Hence 


(5.3) 


i.ryT  ,,2  ^ 2 ~2 

IIVII  - Oi  - CTi 


Combining  (5.1),  (5.2),  and  (5.3)  gives 


Hall2 

2 

Thus  a large  value  of  (o^/5^) 
to  unity. 


(ai^i.)2‘l 

‘ (Oi/o.)2^ 

will  be  reflected  by  the  nearness  of 
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